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Abstract 

In the field of biomedicine magnetic beads are used for drug deliv- 
ery and to treat hyperthermia. Here we propose to use self-organized 
bead structures to isolate circulating tumor cells using lab-on-chip 
technologies. Typically blood flows past microposts functionalized 
with antibodies for circulating tumor cells. Creating these microposts 
with interacting magnetic beads makes it possible to tune the geom- 
etry in size, position and shape. We developed a simulation tool that 
combines micromagnetics and discrete particle dynamics, in order to 
design micropost arrays made of interacting beads. The simulation 
takes into account the viscous drag of the blood flow, magnetostatic 
interactions between the magnetic beads and gradient forces from ex- 
ternal aligned magnets. We developed a particle-particle particle-mesh 
method for effective computation of the magnetic force and torque act- 
ing on the particles. 
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1 Introduction 



1.1 Circulating Tumor Cell (CTC) 

CTCs detach from a tumor and can remain in the blood even after the tumor 
is removed. Their presence increases the chance of net tumors developing. 
It is important to monitor the number of CTCs in the blood but their low 
concentration (a few per /il) when compared to normal blood cells (around 
5 million per fit) makes this difficult. A new and flexible method is required. 

Lab-on-chips with fixed arrays are designed to use the properties of the 
CTCs to filter them. The most common ways are mechanical filters or using 
antibodies. Membranes with slots smaller than tumor cells fulfill the require- 
ments for a mechanical filter [JJ. Normal blood cells, they have similar sizes, 
are more deformable than CTCs and can go through. Another possibility is 
to use cylindrical posts coated with antibodies [2|. If there is a large surface 
area the possibility for a connection with the tumor cells increases. 

Recently a CTC-chip based on guided self-assembly of magnetic beads 
was proposed [3J. Saliba et al. use magnetic traps made by microcontact 
printing in order to create magnetic chains on a regular grid. They demon- 
strate cell capture using beads coated with antibodies. In this paper we 
show that tunable microfluidic chips can be designed by carefully selecting 
the susceptibility and the applied magnetic fields. Using discrete particle 
simulations we compute the distance between particle chains in a microflu- 
idic chip as function of a homogenous magnetic bias field. Such devices might 
be useful to combine immunomagnetic with mechanical filtering of CTCs [I]. 
Maimonis et al. show that the capture efficiency can be increased in a mi- 
crofluidic chip with 2 post arrays of different gap size which enables affinity 
and size capturing. Our proposed chip technology may offer the possibility 
to switch between affinity and size capturing through changing the distance 
between antibody-coated chains by an external field. 

1.2 Content of the paper 

To optimize the design of CTC-chips there are several questions for us. Can 
we use self-organized magnetic beads to create filter-like structures? What 
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types of particles are applicable? Properties to deal with are susceptibility 
and diameter of the beads. To answer this questions we developed a simu- 
lation model with interacting magnetic beads to create particle chains. This 
is described in section 2 of the paper. In addition to the magnetic part we 
developed a blood cell model that flows through the channel with the chain 
barrier. The described methods are used to demonstrate tunability of the 
chain gaps with an external magnetic field. Results are presented in section 
3. 



2 Methods 

2.1 Magnetic particle dynamics 

In order to create a lab-on-chip device consisting of self-organizing micro- 
magnetic beads we want to profit from an external magnetic field and the 
force of the blood flow. Fig. [TJ shows the overview of the proposed microflu- 
idic chip. Viscous blood flows into a pipe filled with soft-magnetic beads. 
Magnetic charge sheets with opposite magnetization orientation provide an 
external magnetic field. 

Under the influence of the magnetic field a magnetic moment m is created 
in every particle. With the moments of two nearby beads and the distance r 
we got a formulation (Eqn. [TJ of the interaction force F t for bead 2 and vice 
versa for bead 1 [5]. 

3fi ( _ , _ _ _ h{rhxr){m 2 r) ^ 

= - — d(mir)m 2 + (m 2 r)mi + mim 2 r - r (1) 

The gradient force F g (Eqn. [2]) on a bead is given by the negative gradient 
of the energy of the magnetic dipole moment m in the field B. 

F g = V{mB) (2) 
In three dimensions it is given by 

(F x \ ( m x dxB x + rriydxBy + m z dxB z \ 
F y J = m x dxB y + rriydyBy + m z dyB z (3) 
F Z J \m x dxB z + rriydzBy + m z dzB z J 
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Figure 1: Formation of particle chains through the interaction of the drag 
force Fd and magnetic forces. In the picture the magnetic forces are cre- 
ated by the magnetic field gradient F g and magnetic interactions Fj. The 
nonuniform magnetic field is created by 2 magnetic charge sheets. 

We assume that in the x-y plane the external field is homogeneous due to 
the distance and size of the charge sheets. For the calculation of the gradient 
force only the z- field derived in y is important. 



(4) 



The blood flow creates a drag force [6] on the particles. At a low velocity 
of 460 [im/s as used at most in CTC-chips [2], the flow is laminar. Turbulent 
flow would appear at Reynolds number (Eqn. [5]) more than unity. 
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During laminar flow the Stokes Law (Eqn. E]) is used to calculate the force 
F d on an object in a fluid, with particle radius r, density p and viscosity r\ of 
the fluid and the relative velocity v of the particle. Krishnamurthy et al. [7] 
show that Stokes drag is a good approximation for the force acting on the 
beads. 

F d = 6nr]rv (6) 

For the above mentioned calculations of the magnetic particle dynamics 
we decided to expand the open-source particle simulator Yade [8]. It pro- 
vides for example gravity force and collision detection using various engines. 
Additionally we developed magnetic and fluidic engines. In the simulations 
we integrate the Newtons equation of motion for magnetic particles under 
the influence of the drag force of the fluid, the gradient force of the external 
magnetic field and the gradient force of the field created by the particles 
themselves. 

For every particle in the pipe we calculate the magnetic field from the 
permanent magnets analytically [9]. We represent the 2 magnets by charge 
sheets as commonly done by magnetic simulation [TO]. The simulation pa- 
rameters are listed in table [TJ 

For blood cell dynamics we implemented a cellular engine in Yade. This 
will be explained in the next section. 



2.2 Blood cell dynamics 

The red blood cell consists of a cytoskeleton covered with a double layer 
membrane [TT]. It has a biconcave form in its native state with a diameter 
of around 8 [im, a surface of 135 /im 2 and a volume of 94 /an 3 . The struc- 
ture is built from 33000 hexagons which in reality would contain different 
proteins such as Spektrin and Aktin. In its loose state there is almost no 
interaction force between these proteins (Fig. [2^). If there is an external 
force acting on the cell the cytoskeleton expands and gets stretched (Fig. 
Hb). The membrane is an almost incompressible layer regarding the plane 
and shear stresses. Capillaries have a smaller diameter than the red blood 
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radius 


KT 6 m 


susceptibility 


0.4 


magnetic saturation 


1 T 


Geometry 




magnetic charge sheets 


10 x 10 mm 


pipe 


10 x 0.05 x 0.02 mm 


distance pipe-sheet 


2 mm 


pole density sheets 


1.6 T 



Table 1: Simulation parameters 



cells, but because of their elastic properties they make their way using the 
form of a projectile. 




Figure 2: Cytoskeleton: a) relaxed, b) stretched 

Discrete element methods can be used to treat red blood cells numeri- 
cally. In the following we will show that the behavior of red blood cells can be 
modeled using a nonlinear mass-spring-system (Fig. [3]). The model consists 
of nearly 600 spherical particles connected by springs arranged orthogonally. 

The force between two particles i and j (Eqn. [7]) depends on the respective 
spring constants, starting distance l and the actual distance d. Yade provides 
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Figure 3: Nonlinear mass-spring model of red blood cell: /c s i=500, k s2 =4:00 } 
fc s3 =300, k s4 =200 (Eqn. [7} 



the damping for all particles in the model. 



F{ Xi , Xj ) =k.*Tjr* (\d\ - l ) * \(\d\ - l )\ (7) 

One possible way of characterizing a red blood cell is using an optical 
trap for extracting its shear modulus [T2] . With an applied force on opposite 
sides of the cell, due to optical tweezers, the diameter is measured. From 
the approximated formula (Eqn. [8]) the modulus decreases form the starting 
diameter D with a function of the force F and the elasticity //. 



To validate this model we computed this cell diameter as function of an 
applied stretch force (Fig. Hj). Different spring constants are used to recreate 
a similar behavior in the stretching test. They decrease from the middle to 
the edge of the red blood cell (Fig. [3]). The modeled cell is accurate enough 
to make tests with the chain barrier. This will be demonstrated in the result 
section. 
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Stretch force [N] 



Figure 4: Stretch test of red blood cell, original [12] vs. model 

3 Results 

3.1 Single particle equilibrium position 

With a steady blood flow through the microfluidic chip the magnetic parti- 
cles are flushed out without the presence of an external magnetic field. To 
find the equilibrium position y of a single bead after applying the field we 
calculate the force balance of the drag force and the gradient force F g 
(Eqn. [9]). Fig. [6] shows this forces and their sum as function of the position. 
At the equilibrium state F^ and F g cancel each other. 



F g (y) = F d =}► y (9) 
The magnetic moment rn(y) of a single particle depends on its volume 
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V, susceptibility Xi which defines the ascending slope of the magnetization 
curve (Fig. [3]), and the magnetic field. This field is the sum of the gradient 
field B g (y) and a homogenous bias field Bf, ias , which we will need to tune the 
device. 



m(y) = V X —(B g (y) + B btas ) (10) 

We have now 2 possible working points for the softmagnetic particles. 
In its saturated state (Fig. [5J1) the magnetic moment can't be further in- 
creased. For better tunability we are using the working point in the initial 
magnetization curve (Fig. [5j2). This will be shown in the next sections. 



ivu-. 




Figure 5: Operation points in the magnetization curve: 1. magnetic satura- 
tion, 2. initial magnetization (tunability) 



The maximum field from the charge sheets is in the center and decreases 
towards the edge (Fig. [6]). Fig. [7] shows the negative drag force and total 
magnetic force depending on the position of the magnetic particles. The re- 
sults show that the final position depends on the strength of the uniform bias 
field and the susceptibility \ of the particles. If there is no crossing point 
the drag force is to high and the particles get washed out of the pipe (Fig. 

with a bias field of less than 0.5T). 

If one of the simulation parameters changes, e.g. the viscosity of the 
blood, the drag force changes and therefore the external magnetic field has to 
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change. Otherwise there probably wouldn't be a crossing point, and therefore 
no equilibrium position. 




Figure 6: Magnetic charge sheets create a gradient field B g . A homogenous 
bias field B g is added to the existing one for better tunability. Force balance 
of fluidic drag force Fa and magnetic gradient force F g leading the particles 
to the equilibrium position y . 



3.2 Chain formation 

The simulation starts with softmagnetic beads of radius lfim and suscepti- 
bility 0.4. They are randomly located in the fluidic channel (Fig. Efet) - Imme- 
diately after applying the magnetic charge sheets the particles self-organize 
to chains according to the field lines (Fig. Hb). This chain formation is then 
shifted to the equilibrium position y because of the force balance explained 
in the section before. Because of chosen fluid velocity of 460 /xm/s this lasts 
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Figure 7: Gradient force F g as function of position with susceptibility a) 0.2 
b) 0.4 c) 0.8. Crossing points with the drag force F d define the equilibrium 
position y of a magnetic particle. 



around 16s (Fig. |8b). 

After this procedure the filtering of the blood could start already. Dif- 
ferent types of CTCs may need different distances between the chains for an 
optimal filtering. In the second phase of the simulation an additional ho- 
mogenous bias field is applied on the chip. Depending on the field the chains 
change their gaps in around 0.13s (Fig. [Hfe). How it works will be explained 
in the next section. 
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a) randomly orientated b) chains created c) equilibrium position 




Figure 8: Timeline of the simulation: Phase 1 - Particle chain creation, Phase 
2 - Changing the gap size 

3.3 Tunable gap size 

Eqn. [TT] shows the magnetic moment rn(y) of a single bead with the exis- 
tance of many particles. When we increase the bias field Bf, ias the magnetic 
moment and therefore the magnetic field created in every particle gets more. 
But only if we are in the working point of the initial magnetization curve 
(Fig. EJ2). The interaction field B t is the sum of all other particle fields. 
This leads to a higher particle interaction force which causes a larger gap. 

Fig. [9] shows the increasing average gap size between the chains. In order 
to quantify this effect we computed the average distance between magnetic 
particle chains as a function of an additional homogenous magnetic field and 
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Figure 9: Average gap size between magnetic particle chains with different 
susceptibilities and additional external field 

susceptibility. 

m(y) = V X —(B g (y) + B Has + B { ) (11) 

With a susceptibility of 0.2 resulting gaps occur only with an external 
field of more than 0.5T. In this case the force caused by the fluid is greater 
than the maximum force of the gradient. So the particles get washed out 
of the pipe. At an intermediate susceptibility of 0.4 the gap size increases 
linearly with the field from 18/im to 24//m. At high susceptibility the total 
field will saturate the particles which leads to similar results after 0.5T. 



3.4 Particle density 

For the simulation result it is essential how dense the particles are filled into 
the pipe. For the particle density Pbead we calculate the number of beads N 
times the volume of one bead Head over the volume of the wrapping box Vb ox 
(Eqn. H2J. 
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(12) 



To evaluate the simulation we compared different starting densities and 
examined the resulting structures. Up to a density of 0.07 we find good 
results with not more than 2 bead thick chains. After that the particles get 
clumpy and can't be used for our application. Fig. [10] demonstrates the 
different results of bead densities. 




Figure 10: Endposition of chains with particle density 0.03, 0.07 and 0.11 



3.5 Cell chain interaction 

Very important for the functionality of the filter is the stability of the chains. 
When a blood cell touches such a chain a certain force is applied on it. And 
it would destroy the structure if it gets to high. To get a value of the force we 
started a simulation of our blood cell that flows directly on a particle chain 
that is fixed in space. Fig. [11] shows the average force per ms over time 
acting from the blood cell to the chain, especially on one magnetic bead of 
the chain. The simulation overview is shown in Fig. [12] The magnitude of 
force on a single magnetic bead in the chain is 1000 times smaller than the 
drag force of the fluid or the gradient force of the magnetic charge sheets. So 
it doesn't influence the position and the stability of the chain structure. 
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Figure 11: Force from the blood cell acting on a single magnetic bead of the 
chain structure. 

4 Summary 

Flexible ways are important to get a high probability of catching cancer 
cells. Micromagnetic beads are established to create chain barriers on de- 
mand. Gap sizes are modified as desired for different CTCs. Size-based 
filtration and the usage of antibodies covered on the beads can take place 
simultaneously. In this work we demostrated simulation tools for the design 
of micromagnetic CTC-chips. 

First step of modelling blood cells was the creation of a nonlinear mass- 
spring-system. For the success of the project it is essential to create models 
of all kind of blood cells, healthy ones and circulating tumor cells. More com- 
plex fluidic behavior is also important for future work. Turbulences could 
occur due to the selfarrangement of the chains in the magnetic field or the 
interaction of magnetic beads and blood cells. 

Results summarized: 
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Figure 12: Red blood cell collides with a single magnetic chain for determing 
the contact force. 

• Starting density of micromagnetic beads shouldn't be more then 0.07 
for the ratio between volume of beads and wrapping box. 

• Softmagnetic beads self-organizes to particle chains due to the external 
aligned magnetic charge sheets within a few fxs. 

• After about 16s the chain arrangement finds its equilibrium position 
with the force balance of fluidic drag force and magnetic forces. 

• Applying an additional homogenous field changes the gap size between 
the chains according to different circulating tumor cells in around 0.13s. 

• A nonlinear mass spring cell model in the flow creates a force on a 
magnetic bead in a chain with a magnitude of 1000 smaller than the 
drag force or the gradient force. So it has no influence on the chain 
structure. 
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